Ramsey interferometry with a spin embedded in a Coulomb chain 
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We show that the statistical properties of a Coulomb crystal can be measured by means of a 
standard interferometric procedure performed on the spin of one ion in the chain. The ion spin, 
constituted by two internal levels of the ion, couples to the crystal modes via spatial displacement 
induced by photon absorption. The loss of contrast in the interferometric signal allows one to 
measure the autocorrelation function of the crystal observables. Close to the critical point, where 
the chain undergoes a second-order phase transition to a zigzag structure, the signal gives the 
behaviour of the correlation function at the critical point. 



I. INTRODUCTION 

The quest for control of quantum dynamics of systems 
with increasing size is one of the present challenges in 
technological applications of quantum mechanics [1]. It 
involves understandingthe transition from the quantum 
to the classical world ^ and is based on the full knowl- 
edge of how the quantum properties scale with the sys- 
tem size, and in particular, of how the thermodynamic 
properties are related to the system microscopic dynam- 
ics [1, 0, @ • Several experiments pursue a bottom- up ap- 
proach, where systems of increasing complexity are built 
by combining together simpler systems, over which one 
has full control In this context, Coulomb crystals of 
ions in Paul and Penning traps constitute a prominent 
system. These crystals are composed by cold ions in a 
confining potential that balances the Coulomb repulsion. 
The ions vibrate around fixed positions in analogy to the 
situation in an ordinary solid, while the interparticle dis- 
tance is usually of the order of several micrometers, con- 
stituting an extremely rarefied type of condensed mat- 
ter 0. Varying the potential permits one to control the 
number of ions, allowing one to explore structures of very 
different sizes, thus offering the opportunity of studying 
the dynamics from few particles to mesoscopic systems. 
Besides, these structures provide promising applications, 
among others, for quantum information processors 0, [1| 
and simulators [i[M[ni,[il- 

Amon g v arious crystalline structures experimentally 
realized jisl \lA [lH . low-dimensional structures, and in 
particular linear chains of ions, are routinely produced 
in ion trap laboratories. Linear chains are obtained in 
anisotropic potentials with large transverse confinement, 
usually linear Paul traps [l^, and can be composed by 
several dozens of ions. They exhibit a mechanical insta- 
bility to a zigzag structure as a function of the ion num- 
ber or of the transverse confinement, which has been first 
characterized experimentally in (TtI . [18| . The dynamics 
and thermodynamics of the linear chain have been ex- 
tensively studied in several theoretical works [H, [13, HH . 
Most recently, it has been demonstrated that in the clas- 
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FIG. 1: Color online. Ramsey interferometry with a chain. 
I. A transverse laser pulse prepares one ion of the chain in a 
superposition of the internal states \g) and |e). The mechan- 
ical effect, associated with the absorption of a laser photon 
and conditioned to the ion being in state |e), displaces the 
ion transversally and excites the modes of the chain. II. The 
chain is let freely evolve for a time t. III. A second laser pulse, 
addressing the same ion of the chain, is sent in the opposite 
direction of the first pulse. The probability that the ion is in 
state \g) is then measured as a function of t. 



sical regime the mechanical instability of the linear chain, 
which leads to the abrupt transition to a zigzag struc- 
ture, is a second order phase transition (2^. [23t. The ef- 
fects of quantum fluctuations on low-dimensional Wigner 
crystals have been discussed in [2^1 , and are expected to 
be in general negligible for Coulomb crystals of atomic 
ions [2^. Most recently, in [2^ parameter regimes have 
been discussed, where quantum effects in structures of 
few ions could be experimentally observed. These theo- 
retical predictions lead to a practical issue, namely how 
to determine experimentally the thermodynamic quanti- 
ties of large crystals. In fact, while the quantum state 
of few trapped ions is experimentally fully accessible by 
quantum tomographic techniques [271 |28|, the applica- 
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tion of such technologies to larger crystals requires di- 
verging experimental times. 

In this paper, we discuss a method for measuring 
the dynamical and statistical properties of an ion chain, 
which is sketched in Fig. [1] This method is based on an 
extension of Ramsey intcrfcromctry with ions [29l |30| . 
which is applied to a dipolar transition of one ion belong- 
ing to the chain. This transition couples to the crystal 
modes via spatial displacement induced by photon ab- 
sorption or emission. By measuring the interferometric 
contrast in a suitable setup, we show that one can di- 
rectly measure the autocorrelation function of the crys- 
tal. In particular, we determine the dependence of the 
interferometric signal on the system size and on the dis- 
tance of the control parameters from the critical point of 
the phase transition to the zigzag configuration. Based 
on the latter dependence, the interferometric signal may 
allow one to measure the critical exponents. 

This work is organized as follows. The salient prop- 
erties of Ramsey interferometry with ions in traps are 
reviewed in Sec.[lTl In Sec. IIIII we introduce the system, 
composed by a Coulomb chain and the spin of one of the 
ions of the crystal. The interferometric signals obtained 
for this system are presented in Sec. IIVI The conclusions 
are discussed in Sec. |Vl In the Appendix we report the 
details of calculations at the basis of the results presented 
in Sec. ITVl 



havior of the loss of contrast in the interferometric sig- 
nal, hence, also permits one to measure some statistical 
properties of the reservoir. This concept has been ap- 
plied for determining quantum properties of cold atomic 
gases [U- 

In this section we review the idea at the basis of the 
setup for ion interferometry, as it was first introduced 
in [30|, and discuss the information the corresponding 
interferometric signal may provide. 



A. Ion interferometry 

We consider a single atom of mass m in a har- 
monic trap, whose electronic states \g) and |e) (both 
metastable) are resonantly coupled by lasers The two 
states |g) and |e) can be, for example, hyperfine states of 
the optical quadrupole transition S1/2 D5/2 in ^"Ca"*" 
whose lifetime is of the order of 1 s Alternatively 
one can couple two metastable states of a hyperfine mul- 
tiplet via a coherent Raman transition using two lasers, 
as in Ref. [l^l, or use two states of a magnetic dipole r.f.- 
transition as proposed in Assuming that only one 
direction of the motion is relevant to the dynamics, the 
Hamiltonian of the system has the form 



int ; 



(1) 



where 



II. RAMSEY INTERFEROMETRY WITH IONS 

Ramsey interferometry is routinely applied in quantum 
optics, nuclear magnetic resonance, atomic and molecu- 
lar physics for high precision measurements [2^ . Most 
recently, it has been considered for measuring the loss 
of coherence in quantum optical systems due to coupling 
with an external environment. Prominent experimental 
applications are in cavity quantum electrodynamics [sTl 
and trapped ions [32]. The idea at the basis of this im- 
plementation can be summarized as follows. Loss of co- 
herence in quantum dynamics is often associated to the 
partition of a large quantum system into subsystems, of 
which only one part S is fully controllable in terms of 
unitary operations and/or of measurements performed, 
while the rest, R, acts as a reservoir. The system S can 
be for instance a two-level transition or a spin, on which 
interferometry is applied. In suitably designed setups, 
one can have spin-dependent coupling to the reservoir R, 
such that the dynamics establish correlations between S 
and R. Such correlations give rise to effects mimicking 
noise and loss of coherence in the dynamics of the ob- 
servables for the spin degrees of freedom Loss of 

coherence can be revealed in the off-diagonal elements 
of the spin density matrix as a function of time, which 
can be measured in the interferometric signal. Their time 
evolution is directly related to the creation of correlations 
with the reservoir R, and thus also indirectly related to 
the statistical properties of the reservoir itself. The be- 
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(2) 



is the Hamiltonian for the internal transition at frequency 
TUe, i.e., the system S, and Hb — hvh^b is the Hamilto- 
nian for the oscillator of the atom center of mass motion 
along the y-direction, i.e., the reservoir R - with h and 
the annihilation and creation operators, respectively, of 
one excitation at energy hv. The term H^t describes the 
interaction of the two-level system with a laser pulse at 
frequency ujl and wave vector along the y-direction, 
and reads 



hn{t) 



-HuJLt-kLV) 



H.c. 



(3) 



where a = \g) (e|, its adjoint, and il{t) is the real- 
valued Rabi frequency. From now on we take ujl = ^0 
and consider the dynamics in the reference frame rotat- 
ing at the laser frequency. In this frame, the phase of the 
field depends solely on the center-of-mass position y of 
the ion. By absorbing/emitting a photon, the atom expe- 
riences a mechanical displacement due to the field phase 
gradient over the center-of-mass wave packet. This me- 
chanical action is described in Eq. ^ by the operator 
exp(iA;iy), which, for the harmonic motion, is the dis- 
placement operator 

D{a) = exp(a6^ - a*b) , (4) 

where we used y = ao(fe^ + b) and a — i/cLOo with gq = 
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We now assume that at time t = the ion is prepared 
in the internal and oscillator ground state |5'(0)) = \g, 0). 
In the time interval [0, tpuise] a square pulse at constant 
intensity Q{t) = Qq is applied, such that Slo'^puise = 7r/4 
(thereby implementing a so-called tt/2 pulse). In the 
regime in which the pulse intensity is sufhciently large, 
such that i^ipuiso ^ 1, the free evolution of the oscillator 
can be neglected during the pulse and the ion is prepared 
in the superposition 



|*(<pulsc)> 



V2 



(|5,0) + |e,a)) , 



(5) 



where |a) = D{a) |0) is a coherent state. 

The system is then let freely evolve for a time t while 
a unitary operation is applied, which introduces a spin- 
dependent phase (/) according to the procedure 
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(6) 
(7) 



and that the interferometric procedure aims at detect- 
ing. At time t + ipuisc, assuming a purely Hamiltonian 
dynamics, the state of the system reads 



1 



|«'(i + ipuisc)> = ^(l5-0) 



e, ae 



*» • (8) 



A squared laser pulse is then switched on for a duration 
Sipuise, thereby implementing a —tt/2 pulse on the spin. 

In the absence of coupling with the motional degrees 
of freedom, the probability to find the atom in the initial 
state after this pulse would solely depend on the phase (j). 
However, due to the coupling with the external degrees 
of freedom, there is a dephasing which arises from the 
free evolution of the coherent state \a). As a result, the 
probability to measure the atom in the ground state after 
the last pulse takes the explicit form 



Re{e''^5(t)}] 



(9) 



with 



Sit) 



(10) 



and < \S{t)\ < 1. For fixed time evolution t, the prob- 
ability Vg as a function of is hence a periodic signal 
exhibiting fringes with visibility 



V^\Sit)\. 



(11) 



In [32| a version of this interferometer has been used in 
order to measure decoherence of the ion quantum mo- 
tion, due to the coupling to a classical noise applied at 
the electrodes of a Paul trap. The interferometric sig- 
nal hence gave a measure of how coherence decays as a 
function of time for different statistical properties of the 
applied noise. 



B. Interferometric signal and reservoir properties 

Let us now discuss the physical meaning of the quan- 
tity S{t) in Eq. From an interferometric point of 
view, we notice that loss of visibility is found at times t 
that are different from an integer multiple of the oscilla- 
tion period of the ion. At these instants, then, the — 7r/2 
pulse does not bring the spin back to the initial state, 
and there are residual correlations between ion motion 
and spin which gives rise to a diminution of the contrast. 
In other words, a which-way information is left in the 
reservoir R after it interacted with the system. This con- 
nection has been elegantly elucidated by Englert in |3Q] . 
where a physical quantity, the distinguishability V, was 
introduced in order to quantify the amount of which-way 
information left by means of a unitary evolution cou- 
pling spin and reservoir. Englert showed that visibility 
and distinguishability are related by the inequality [36| 



(12) 



where the equality holds when the reservoir is initially in 
a pure state. 

We also notice that the visibility of the signal in 
Eq. (fTTj) can be rewritten as 



V(i)= (0|e-''=^2'(*)e^'=^^(°) |0) =cxp 
where 



(13) 



(14) 



is the autocorrelation function of the reservoir, and the 
mean value (.) is taken over the initial state of R. The 
visibility is hence a direct measure of the autocorrelation 
function of the reservoir R. This result holds in general, 
when the state of the reservoir R is Gaussian. Designing 
a different setup will give access to different correlation 
functions of the reservoir. By implementing a mechani- 
cal excitation as in for instance, one can measure the 
anti-symmetric part of the correlation function. Alter- 
natively, by applying the pulses to two different ions one 
can determine the correlation length of the crystal. In 
general, the Ramsey signal will allow one to characterize 
the thermodynamic properties of a reservoir R in thermal 
equilibrium. 



C. Discussion 

The visibility of the interference signal is reminiscent 
of the Loschmidt echo as defined in [37[. There, a quan- 
tity is evaluated which is the overlap between two states 
evolved from the same wavefunction under the influence 
of two different Hamiltonians, and which is analogous to 
S{t) in Eq. Such an overlap is directly related to 

the Loschmidt echo, which is frequently used to describe 
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loss of coherence as a consequence of the interaction be- 
tween a system and its environment. The Loschmidt echo 
has been extensively analyzed for one-dimensional spin 
chains, where the system is one spin and the decohering 
environment is the rest of the chain, constituting a spin 
bath and characterized by Ising or Heisenberg nearest- 
neighbor interactions [ssl- Differing from that model, 
the present work is concerned with an environment char- 
acterized by long-range interactions, which emerge from 
the Coulomb repulsion between the ions of the chain. 
However, we shall show that the Loschmidt echo in jsl] 
and the visibility of the Ramsey signal discussed in this 
work exhibit important similarities when the correspond- 
ing environment is close to a critical point, marking a 
second-order phase transition. 



III. A SPIN COUPLED TO A COULOMB 
CHAIN 

The theoretical model, which we analyze in this work, 
is a Coulomb chain coupled to the two-level electronic 
transition of one of the ions composing the chain, which 
we will denote as the spin. Interaction between the spin 
and the chain occurs via the mechanical effects of light: 
The ion displacement, due to photon absorption or emis- 
sion, induces a mechanical excitation of the chain modes 
via the Coulomb and trap forces binding the crystal. The 
Hamiltonian of the system reads 



(N) 



int : 



(15) 



where i?at is the spin Hamiltonian given in Eq. 
describes the dynamics of the chain, and iJjnt gives the 
interaction of the spin with the laser pulse and is defined 
in Eq. ([3]) . In the first part of this section we review the 
basic properties of the ion chain. In the second part we 
characterize the mechanical effects of light on the crystal 
due to the coupling with the laser pulse, described by the 
term Hi^t. 



with periodic boundary conditions, where the ions ex- 
perience a mutual repulsive interaction along the linear 
axis, which we identify with the a;-axis, while the trans- 
verse direction is confined by a steep harmonic potential 
of frequency vt [soj . We assume that one ion is pinned at 
the origin, so that the equilibrium positions of the other 
ions arc rj"' = (x^p ,0, 0) with = ja, j = 0, 1, . . . , iV, 
and a is the interparticle distance, a ~ Lchmn/N. More- 
over, we assume that the ions are sufficiently cold such 
that their vibrations around the equilibrium points are 



harmonic. We denote by qj 



„(") 



Uj, and Zj the 



axial and transverse displacements of the ion j from its 
equilibrium position. Their dynamics are then described 
by the Hamiltonian 

H^'''> = Vo + H.,,+Hy + H,, 

where Vq denotes the energy of the classical ground state 
and Hi (I = x, y, z) give the harmonic motion about the 
classical equilibrium positions. In the normal modes de- 
composition, after quantizing the linear excitations, these 
terms take the form [2l| 



k <J=± 



(16) 



where k labels the mode wave vector, which takes the 
values k = 2'Kn/Na with n = 0, . . . ,N/2, while cr = ± 
denotes the mode parity (see also Appendix A) . The op- 
erators bi{k, cr) and bl{k, a) annihilate and create, respec- 
tively, a quantum of energy hu!i{k) of the corresponding 
oscillator, whereby |4q| 



= 1 



(17) 



JV/2 



u.,{kf = uj.{kf = v\ - ^u^l Y: 4 ^ - (18) 



with 



A. The Coulomb chain 

The Coulomb chain is composed by N ions, each hav- 
ing the same mass m and charge Q, and confined inside 
a highly anisotropic trap, which stabilizes the crystalline 
order along a line'^. For convenience, in the rest of this 
paper we will consider a ring trap of very large radius. 
In this regime the system is equivalently described by a 
linear distribution of charges on a line of length Lchain 



^ Experimentally, ion chains are usually realized in linear Paul 
traps. Some experiments reported ion chains in ring traps, as for 
instance in llTll . Here, crystallization is achieved combining laser 
cooling, to the regime in which the ions thermal energy is much 
smaller than the interaction energy, with some sort of pinning. 



(19) 



a characteristic frequency scaling the spectral frequen- 
cies of the crystal. We notice that the transverse trap 
frequency Vt is the largest eigenfrequency of the trans- 
verse modes. 

The displacements qj , yj , Zj of the ion at position rj°'' 
are related to the normal modes by the orthogonal matrix 
7^, such that [III, [U 



Vs = 




2mujy{k) 



7^,,fc. [6.(fc,fT) + 6t(fc,^)] (20) 



■Rj^ka [hy{k,a)+hl{k,a)\ (21) 
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and analogously for Zj. The specific form of the orthog- 
onal matrix TZ is given in Appendix \^ 

The mechanical stability of the ion chain is warranted 
provided that the transverse confinement satisfies the re- 



the zigzag depends on the value of vt and it is reported 
in Appendix VK\ 

In this work, we will also be interested in studying the 
interferometric signal as a function of the distance of the 

fc) (c) ... . . . (c) 

lation > , where Vf is a critical frequency which transverse frequency from the critical value . We will 



depends on ujq and on the density of ions 1/a. For ions 
on a ring one finds [12, 



(22) 



At i^t 



the chain undergoes a second-order phase 

(c) 

transition to a planar structure, such that for z^t < 
the ions organize in a zigzag structure across the a:- 
axis. Assuming that the zigzag structure is localized 
on the X — y plane, the new equilibrium positions are 



0) 



{x'^p ,yf\Q) where x 



and 



,(0) - 



(0) 
3 



ja, = (-1)^6/2, 

0. The transverse displacement b depends 
on the linear density of ions a and on the value of vt, 
as shown in Ref. [23|. In the zigzag configuration, the 
Hamiltonian describing the harmonic oscillations around 
the equilibrium positions takes the form 



(N) 



TJZZ , TJZ 



where V™ is the energy of the new classical ground state. 
The Hamiltonian is the harmonic term for the dis- 
placements qj and Wj = yj — y^ around the new equilib- 
rium positions, which are now coupled. The dispersion 
relation is now defined in the Brillouin zone [0,7r/2a], 
such that the wave vector k takes the values k = 2iTn/Na 
with n = 0, . . . , N/A. The spectrum of the excitations 
along X and y presents four branches, which we label by 
/3 = 1,2,3,4 (a detailed discussion of the form of the 
spectrum can be found in 23]). In the normal modes 
basis Hamiltonian H^^., reads 



HZ 



k (T=±/3=l 



4i{k,cr)c,3{k,(j) 



(23) 



where a is the mode parity. Finally, the operators 
Cf3{k,a) and c|j(fc,CT) annihilate and create, respectively, 
an excitation of the mode of the branch f3 with wave vec- 
tor k, frequency to/sik), and parity a. Hamiltonian 
is the harmonic term for the displacement Zj and is not 
explicitly reported, as it will not be relevant for the rest 
of this paper. 

The displacements qj, Wj are related to the normal 
modes by the orthogonal matrix TZ™ according to the 
relations 



then denote this quantity by the parameter 



(25) 



which has the dimension of an angular frequency. In par- 
ticular, when A > the linear chain is stable, while for 
A < the system is found in the zigzag configuration^. 



B. Ramsey interferometry with the Coulomb chain 

Ramsey interferometry with the Coulomb chain is im- 
plemented by means of a straightforward extension of the 
procedure described in Sec. Ill Al for a single ion. The ba- 
sic setup discussed in this paper is sketched in Fig. [TJ 
We assume that the ion at position ri in the chain is 
selectively addressed by the laser pulses, which propa- 
gate perpendicularly to the chain, and restrict the Hilbert 
space of the electronic degrees of freedom to the two level 
of its internal transition. By absorbing and emitting a 
photon during the pulse, the ion undergoes a transverse 
displacement, which excites all modes of the chain cou- 
pling with the ion position. As in Sec. Ill Al we assume 
that the duration tpuiso of the laser pulses is sufficiently 
short, such that the free evolution of the chain during 
the pulse is negligible. This corresponds to imposing the 
relation Wmaxipuiso ^ 1, where Wmax is the largest eigen- 
frequency of the crystal modes that are involved in the 
dynamics. In this paper we consider laser pulses propa- 
gating along the y-direction, inducing hence mechanical 
displacement along y. Hence, when the ions form a lin- 
ear chain Wmax = i^t and the condition to be fulfilled is 

The ion displacement due to absorption or emission of 
a photon is described by the operator 

exp(ifcL2/i) = JJexp [i7y(fc,(T)(&j^(fc,(T) -t- 6|^(fc,(T))] 

k.a 

= llDkAc^ka), (26) 

which has been here conveniently expressed in terms of 
the normal modes of the linear chain using Eq. (pij) . Here 



rj{k,a) = kLi 



2mujy{k) 



k(T 



(27) 



k cr=± 13 



2mujf{k) 



T^TkBa \c0{k,a)+cl{k,a) 



(24) 

where g — (gi, wi, . . . , qjv, fi'Ar) and the index j = 
1, . . . , 2N. The specific form of the orthogonal matrix for 



The structure is found in a zigzag configuration for values of ut , 

(c) (c) 

such that ' > Ut > i^l 



where ' is another critical value 



at which the system undergoes a phase transition to a multiple 
chain structure [13, SHI- 
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is the Lamb-Dicke parameter for the mode tOy{k) scahng 
the displacement of the mode by photon absorption or 
emission [ill, and Dk^a{aka) denotes the displace- 
ment operator for the oscillator at frequency LOy{k) and 
parity a with amplitude 

aka = cr) . (28) 

For later convenience we denote by 

\{ak,a}) =\{DkAakam = ^\aka{t)) (29) 

k,a k.a 

the quantum state of the linear chain, obtained by apply- 
ing the displacement operator (1^^ to the ground state 
of the chain |0). 

Let us now consider some useful relations. The chain 
will be in the Lamb-Dicke regime, i.e., the ion dis- 
placement by photon recoil perturbs weakly the me- 
chanical state of the chain, provided that the relation 
\ri{k, a)\y^2nka + 1 ^ 1 is satisfied for all modes, where 
nka = {b\(k,a)hy{k,a)) is the mean occupation of the 
mode with frequency Wj^(fc). Correspondingly, \aka\ ^ 1- 

In general, in the rest of this work we will assume rjy <^ 
1 , where rjy — rjo/ y/N is the Lamb-Dicke parameter of the 
transverse bulk mode at frequency i^t , and 



2mvt 



(30) 



while we will not make any particular assumption for the 
Lamb-Dicke parameter of the other modes. In particular, 
close to the mechanical instability to the zigzag config- 
uration, the Lamb-Dicke parameter of the modes whose 
frequencies vanish, ujy{k) 0, can become very large. 
Hence, in such a regime, the photon recoil can perturb 
significantly the state of the chain. 

Finally, when studying the interferometric signal we 
will compare curves obtained at different values of vt- In 
this case we will make use of the quantity 



ic) 



2mvl 



(31) 



which is related to the Lamb-Dicke parameter 770 by the 

equation rj^'^^ — rjQ \J vt/v['^^ . We conclude this section 
by giving some typical values of the parameters that are 
relevant in this work. In the experiment of Ref. [l3], 
ion chains of dozens of ^^Mg"*" ions were realized, with 
interparticle distance of the order of a '--^ 33 /zm. This 
case corresponds to the value of the typical frequency 
Wo ~ 27r X 64 kHz, while the single particle Lamb-Dicke 
parameter, close to the critical point, is 770 ~ 0.62. 



IV. MEASURING THE STATISTICAL 
PROPERTIES OF THE CRYSTAL 

We now analyze the signal obtained by performing 
Ramsey interferometry with the spin of the ion at the 




(n/co„ 



FIG. 2: a) Visibility V as a function of the time elapsed 
between two Ramsey pulses, in units of 1/ljo, for A'^ = 100 
ions in a linear chain. The parameters are ut = z^^'^' -I- A, 
with A — lQ~^Loo, and r;''^' = 0.25. b) Fourier spectrum J-{lj) 
of the visibility V(t), obtained by computing numerically the 
discrete Fourier transform of V{t) in a finite time interval 
t e [-Tf/2; Tf/2] with Tp = lO^/oJo- The interval is sampled 
with a time slicing dt = Tf/tis with Us — 10^. We checked 
convergence of the spectrum amplitude JT by increasing Tf 
and Us. The spectrum is normalized such that J^{uj = 0) = 
1. The dashed lines indicate the minimum and maximum 
frequencies of the transverse spectrum (tJmin = y/2/^Vt + A2 

and CiJmax = Vt). 



position ri of the chain. We assume that the spin is ini- 
tially in the electronic state \g) and that the chain has 
been prepared in the ground state by some cooling pro- 
cedure [3a|- Hence, the initial state in the Hilbert space 
of the modes of the linear chain and of the spin transition 
reads 



l*(0)) = |g,0) 



(32) 



After the sequence of Ramsey pulses, chain and spin are 
in the entangled state 



|*(t)) = 1^,(0) \g)+i\^M)\e) 



(33) 



where \ipg) and |(/?e) are (not normalized) states of the 
crystal and read 

- \ m+^-'^''^'^\{o^Um] , (34) 
We{t)) - i[|{«fc4>-e'^|{«fc.(i)})]. (35) 
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FIG. 3: Same as in Fig. [2] except for vt = ^'i + ^ with 



Here, aka{t) = afco-e"''^''^'^^* and (j) is the phase due to the 
appUed phase shift. The visibihty in the Ramsey signal 
is 

V = |({«fe.}|{afca(i)})| = exp[-A(t)] , (36) 
where, for a hnear chain, the exponent reads 

A(t) = 2^K.Psin^^^^. (37) 

k{7 

An analogous expression can be obtained for the zigzag 
structure. The quantity A{t) is directly related to the 
autocorrelation function in Eq. (|14p . 

A{t) = ^g{t), (38) 

where this relation is valid both for the linear chain and 
for the zigzag structure. Moreover, it also holds when 
the crystal is in thermal equilibrium at temperature T. 
In this case the function A takes the form 

A^{t) = 2;^coth^|a..Psin^ ^ , (39) 

k{7 

where fcs is the Boltzmann constant. In the rest of the 
paper we will focus on the case T = for which the 
second-order classical phase transition linear-zigzag at 
Vt = vl occurs. 



Figure [2^) displays the visibility V for the linear chain 
as a function of the elapsed time t for N = 100 ions and 
7](^^ = 0.25. The chain is sufficiently far away from the 
critical point (A = 10~^ijJo), such that all modes arc in 
the Lamb-Dicke regime and the mechanical effect of the 
photon perturbs weakly the chain stability. Nevertheless, 
we observe that the visibility exhibits initially a decay, 
reaching then values below unity, about which it oscil- 
lates asynchronously for t > Q. This behaviour is a clear 
consequence of the time-dependent form of the signal, 
as it can be seen expanding Eq. ([55]) in powers of the 
exponent A{t), defined in Eq. ([57)1 : Hence, the Fourier 
spectrum of the visibility contains in principle all possi- 
ble sums of the eigenfrequencies (^y{k). The dispersion 
relation shows that the eigenfrequencies are incommen- 
surate, such that the signal in Eq. ([57)1 is not periodic (al- 
though for finite systems it may exhibit partial revivals 
in the excitations, as we will discuss in Sec. IIVBI) . Fol- 
lowing the interferometric interpretation by Engiert [36j , 
the mechanical effect of light leaves a which-way infor- 
mation in the chain excitations, which is only partially 
erased by the second Ramsey pulse. 

When the structure is in the Lamb-Dicke regime, the 
visibility signal can be put in direct connection with the 
autocorrelation function g{t), since V{t) ~ 1 — A{t). The 
Fourier spectrum of the visibility, in this case, is thus 
just composed by the frequencies of the normal modes. 
Figure [2b) displays the Fourier spectrum of the visibility 
in Fig. [5^), obtained by performing the discrete Fourier 
transform of the signal V(^). The sharp peaks, giving the 
main components of the Fourier spectrum, correspond 
here to the frequencies LOy{k) of the normal modes. 

Let us now discuss the signal when the parameters of 
the linear chain are close to the critical point, where it 
undergoes a transition to a zigzag configuration. Fig- 
ure [3] displays the visibility and its Fourier spectrum for 
A = lO^^cjo- The visibility signal decays to smaller val- 
ues as a function of the time elapsed between the two 
Ramsey pulses, thereby exhibiting fast oscillations about 
the decaying mean value. The Fourier spectrum shows 
that the low frequency components of the transverse nor- 
mal modes mostly contribute to the signal. In particular, 
the lowest frequency mode, which corresponds to the soft 
mode driving the instability to the zigzag [1^, has the 
largest Fourier amplitude. In fact, even if the chain is 
still in the Lamb-Dicke regime, the Lamb-Dicke param- 
eter of the soft mode is the largest, thus this mode is 
mostly excited by the ion displacement due to photon 
recoil. In this case, hence, the visibility allows one to ac- 
cess the behaviour of the autocorrelation function close 
to the critical point. 



A. Short times: quadratic decay of the visibihty 

We now focus on the behaviour of the visibility for 
short elapsed times t. In this regime the autocorrelation 
function g(t), and correspondingly the logarithm of the 
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FIG. 4: Color online. Function A{t) as a function of the time 
elapsed between two Ramsey pulses (in units of I/loq). The 
solid line corresponds to Eq. p7|) and the symbols correspond 
to the short-time expansion in Eq. (|4Up . The curve has been 
evaluated for r;''^' — 0.05 in a linear chain of = 10'^ ions. 

visibility A{t) in Eq. ([57)1 . can be expanded in powers of 
u!y{k)t ^ 1. At lowest order in the expansion one finds 

A~Tt^ , (40) 
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FIG. 5: Color online. Coefficient F (in units of ajg) as a func- 
tion oi /\ = vt — vj:'^^ (in units of loo). The solid line corre- 
sponds to the coefficient F, given in Eq. (|4ip and the symbols 
correspond to the value Fflt, obtained by fitting the function 
A{t) in Eq. (|37p with the function Fflt^^- The parameters are 
the same as in Fig. |4] Inset: Derivative of F with respect to 
A (in units of ujo) as a. function of A (in units of luq and log- 
arithmic scale) for A > (symbols). An analytical analysis 
of the asymptotic behaviour, as well as a numerical fit of the 
curve give a logarithmic dependence of the derivative on A 
for A ^ 0+ (line). 



where 

r = i^K.pa;,(fc)^ (41) 

k,(7 

and the condition Wmax^ ^ 1 must hold. The correlation 
function hence goes quadratically with time, as shown in 
Fig. m and correspondingly the visibility signal decays 
with a Gaussian-type behaviour. In particular, for the 
linear chain the coefficient F can be rewritten as 

and is hence proportional to the mean value of the fre- 
quency of the transverse excitations. The detailed deriva- 
tion of this expression is reported in Appendix [B] Equa- 
tion (|42p has been obtained for ions in a ring trap with 
large radius. It approximates well the result for a linear 
chain in a linear Paul trap when N ^ 1. 

We now analyze the dependence of the coefficient F 

(c) 

on vt for values close to i/( , thereby getting further in- 
sight on the behaviour of the correlation function close 
to the critical point. Figure [5] displays F as a function of 
A = Vt — vf^ for values close to the value A = 0. One 
clearly observes a minimum, with a cusp-like behaviour, 
at the critical value vl corresponding to A = 0, show- 
ing that decay of the visibility signal as a function of the 
elapsed time is slowed down close to the critical point. 
The discontinuity of the derivative of F with respect to 
A can be attributed to the structural phase transition 
that the chain undergoes at A = 0. In particular, an 
analytical study shows that for A > (on the side of the 
linear chain) the derivative dF/dA ~ InA as A — > 0+. 



This result, whose derivation is reported in Appendix iBl 
is confirmed by a numerical fit, and is shown in the inset 
of Fig. U 

As anticipated in Sec. Ill C[ analogous features have 
been found in the Loschmidt echo of a spin coupled to 
a quantum spin bath, when studying the short-time be- 
haviour of the echo signal as a function of the control 
parameter close to the critical point of the spin chain 
(see in particular Figs. 2 and 3 of Ref. [sst). 



B. Long times: revivals and asymptotic value of 
the visibility 

We finally investigate the behaviour of the visibility 
V(t) when long times t are elapsed between the Ram- 
sey pulses, focusing on values of the trap frequency vt 
close to the critical point. We consider time scales, such 
that the elapsed time satisfies the relation Wmin^ ^ 1- 
While the asynchronous oscillations of the modes of the 
chain, excited by the photon recoil, lead to a decay of the 
visibility as a function of t, in finite systems we expect 
to observe characteristic elapsed times t* where a quasi- 
synchronization of the motion of the chain will occur, 
and, correspondingly, a sort of revival in the behaviour 
of the visibility signal is observed. The time t* can be 
estimated by considering the propagation speed at which 
the mechanical excitation propagates through the crystal 
and returns to the initial position. This corresponds to 
the approximated formula fll| 
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FIG. 6: Color online. Behavior of the visibility signal for 
long times. The exact result for the visibility V (solid) ob- 
tained using Eqs. (|36p and (|37|l and the approximate function 
exp[—Aoo — B{t)] (symbols) (see Eqs. (|47)) and (|48l)) as a func- 
tion of time (in units of 1/ujo) for A = 10~^ti;o (linear chain) 
with iV = 10^ and ry^"' = 0.25. The vertical dashed line cor- 
responds to an estimate of the revival time according to the 
formula in Eq. (|43|) . which gives t* — 1229/ii^o. Inset: the 
same plot in a smaller region for a closer comparison between 
V and the approximated expression exp[—Aoo — B{t)]. 



where Vmax = maxfc dujk / dk is the maximum group ve- 
locity in the interval k e [0;7r/a]. Figure [6] displays the 
visibility signal as a function of the elapsed time for a 
linear chain composed by TV = 10'^ ions and close to 
the zigzag instabihty, A = lO^'^o^o. For ■q^'^^ = 0.25 the 
group velocity is maximum at fc = 2.64/a, taking the 
value Vmax = O.Slawo- Using this result, the estimated 
time from Eq. is t* — 1229/wo which approximately 
coincides with the time at which the visibility exhibits 
sudden oscillations of larger amplitude. 

We now focus on the regime, in which the system is suf- 
ficiently large, ^ 1, and derive the asymptotic form of 
the visibility signal for times t —^ oo, obtained by aver- 
aging over time intervals T ^ t*. The asymptotic form 
of V(t) of can be singled by observing that the exponent 
A{t) in Eq. ([55]) can be rewritten as 

A^A^^ + B{t) (44) 

where 

Aoo Mit)) + [vim] = E ' (45) 

k(7 

B{t) = fciRe(yi(i)2/i(0)) =El"fc-l''=os'^2'(^)*(46) 

kcr 

Using Eq. ((38)) we see that A^o is proportional to the 
variance of the ion displacement {yHt)) = (2/1 (0)) around 
its equilibrium position, while B(t) is proportional to the 
real part of the crystal's correlation function {y{t)y{0)). 

For TV ^ 1, and when the linear chain is close to the 
instabihty, A ^ 0+, Eqs. (gg) and (glD can be approxi- 



0.02 




FIG. 7: The asymptotic value of the correlation function, as 
given by Aoo in Eq. (|45p . as a function of A when the linear 
chain is close to the critical point, where it makes a transition 
to the zigzag structure. The symbols correspond to the curve 
given by Eq. (|45|l and the solid line is the curve obtained 
from the analytical prediction in Eq. (|47|) . The parameters 
are iV = 10^ and r;'-'''' ~ 0.05. 

mated by the quantities 

^ ioo = -4^1nA + c, (47) 
Bit) ^ m = -it^NoiSt), (48) 

where c is a constant, Nq is the Bessel function of the 
second kind (see Appendix [B|) and we introduced the 
parameters h = Vln 2 ujoa and 

5 = \[v[^\ (49) 

which is related to the parameter A by the relation 

5 = \J /S.{2v['^^ + A). In Fig. [6] we compare the exact 
result of the visibility with the approximated function 
V{t) — exp[— — B{t)]. For long times, such that 
t ^ 1/(5, the function B{t), and as a result the two- 
time correlation function of the crystal, decays as l/^/Si 
and the visibility approaches the asymptotic value given 
by V^oo = exp[— Aoo]. Figure [7] displays Aoo as a function 
of A, showing that Aqo (and therefore the spatial width 
of the ion wavepacket) diverges logarithmically with ut 
close to the critical point. As a consequence, the asymp- 
totic value of the visibility Voo becomes smaller when the 
ion chain is close to the mechanical instability, where it 
undergoes a second-order phase transition to the zigzag 
structure. This behavior of the visibility is qualitatively 
similar to the behavior of the asymptotic Loschmidt echo 
£00 for a spin interacting with an Ising chain close to the 
critical magnetic field [381] . 

V. CONCLUSIONS 

The aim of this work was to understand the dynami- 
cal properties of a quasi-one-dimensional many-body sys- 
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tern externally driven and probed via one of its com- 
ponents, in a setup that is both experimentally feasible 
and not previously studied in this context - a chain of 
trapped ions. The ions form a Coulomb crystal under- 
going a second-order phase transition between a linear 
and a zig-zag configuration, depending on the transverse 
trapping strength. Building on interferometric methods 
with trapped ions ^3^], we proposed a Ramsey experi- 
ment to measure the autocorrelation function of an ion 
crystal by only exciting and detecting the internal state 
of one of its constituent ions. 

The central quantity to this study is the visibility of 
the fringes of this specific implementation of the Ramsey 
interferometer. We show that it provides information of 
the dynamical and statistical properties of the crystal. 
In particular, for the specific realization here proposed 
the visibility allows one to measure the autocorrelation 
function of the crystal, and the characteristic decay time, 
thereby accessing its properties close to the critical point. 

It is interesting to remark that several features of the 
Ramsey signal, close to the critical point where the linear 
chain undergoes a transition to a zig-zag configuration, 
share several analogies with the Loschmidt echo evalu- 
ated for studying the decoherence of a spi n coupled to 
quantum spin-baths close to criticality [H, IH, |3, \^ . 
In fact, the crystal can be regarded as a reservoir with 
respect to the single spin constituted by the two internal 
states of the probed ion. In this sense our system can 
be seen as a model for decoherence induced by a bath at 
T = 0. On the other hand, a distinguishing feature of our 
proposal is the presence of long-range interactions, typ- 
ically absent in (nearest-neighbor interacting) spin sys- 
tems. 

We conclude by observing that the ability to drive 
and probe the whole crystal by addressing just one ion 
opens interesting questions like the possibility of using 
the methods developed here to identify control proce- 
dures, of spin-echo type, to maintain coherence in the sys- 
tem, as well as to measure other quantities, such as, e.g., 
coherence length and quantum fluctuations at a quantum 
critical point. This points to the need for a thorough un- 
derstanding of quantum phase transitions in ion traps, 
as initiated in [26] . The adaptation and further develop- 
ment of the ideas proposed here to the latter context will 
be the subject of our future work. 
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APPENDIX A 

In this appendix we provide more details on the eigen- 
modes of the linear and zigzag structures. A more ex- 
tensive treatment can be found in (23j . For a linear ion 
chain in a harmonic trap, the reader is referred to [2l| . 

Linear chain. For a given frequency ujx{k), ojy{k) there 
are in general two modes of opposite parity with respect 
to the transformation qj — > q-j, yj — > JJ-j, even a — + 
and odd a — —. The frequencies at wave vectors k ^ 
and k — n/a are instead not degenerate. The correspond- 
ing modes are the bulk and the zigzag mode, respectively, 
which have definite parity, the first being even and the 
latter odd. For the transverse excitations, the bulk mode 
is at frequency Vt while the zigzag mode is at frequency 
ujy{k = IT /a). For the modes at wave vector < fc < ir/a, 
the elements of the orthogonal matrix, connecting the 
displacements of the ions from the equilibrium positions 
with the eigenmodes amplitudes, read 



cos{jka) , 
sin(jfca) , 



(Al) 
(A2) 



while the elements connecting the displacements from 
the equilibrium positions with the bulk and zigzag eigen- 
modes take the form 



^.,0 = 

= (-i)^' 



(A3) 
(A4) 



Zigzag configuration. In the zigzag structure the dis- 
persion relation is now defined in the Brillouin zone 
[0, 7r/2a], such that the wave vector k takes the values 
k = 2'Kn/Na with n = 0, . . . , iV/4. As shown in ^ the 
spectrum of excitations consists of four branches that we 
label j3 = 1,2,3,4. For a given frequency u!^^{k) for 
< k < 7r/2a there are two modes with opposite parities 
a ~ ± which are a combination of displacements in the 
longitudinal and transverse directions. The correspon- 
dent elements of the matrix TZ^^ read 



'R-2]-i,kf3+ = upcos{kai), 
'^2j,k0+ = W/3 sin[(7r - fca)j] 



7^: 



■2j-l,kf3- 



7^: 



2j,kl3 



— u/3 sin[A;aj] , 
- = -vp cos[(7r - ka)j] , 



(A5) 
(A6) 
(A7) 
(A8) 



with j = 1,2,...,N labeling the ion displacements 
P2j"_i — qj and g2j = Wj\ and k = k for (3 = 1,2, 
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k = TT / a — k ior f3 = 3, i. The explicit form of the coeffi- 
cients and V(3 can be found in (23|] (check Eq. (20)). 

For fc = the normal modes are the bulk modes in the 
X and y directions. These modes correspond to the zigzag 
structure oscillating rigidly in the x or y directions. The 
matrix elements for the bulk mode in the x direction are: 



2j-l,0x 



0; 



(A9) 



while for the bulk mode in the y direction they are 



7^: 



2j-1,0jj 



0; 7^: 



(AlO) 



where 5 



,(c)2 



2v^t'''>A and h 



'In 2LUoa. 



Taking the continuum limit in Eq. (jB2[) we find 



Ti/a 



Z TT 



1 2 a 1 

4 TT h 



S arcsinh 



dkujy{k) 



\/d^ + K^q^dq + c 



hqmi 



(B4) 



For k = 7r/2a the normal modes are the two zigzag 
modes along the x and y direction. The zigzag mode in 
the X direction is the mode where neighboring ions oscil- 
late around the equilibrium positions along the x axis and 
with opposite phase. The corresponding matrix elements 
are given by: 



7^: 



2i-l„ 



7^: 



2j, 



0. 



(All) 



Analogously, for the zigzag mode in the y direction, the 
matrix elements are: 



— n- 7?^^ 



(A12) 



where ^max is a cutoff and c is a constant which depends 
on the cutoff. The derivative of F with respect to A, for 
A 0+ is obtained from Eq. (|B4[) and reads 



dr 1 9 Vf'^^'a , , 
— r = — Vni't— In A 

dA 4 '° * TTh 



Ci 



(B5) 



max and and we neglected terms 
0. 



where ci depends on q, 
which goes to zero as A 

Long time behaviour. We now study the behaviour of 
the function A{t) for long times. The quantity Aaa can 
be cast in a simplified form for the linear chain. Using 
the same technique which leads to Eq. (|B2p . one finds 



APPENDIX B 

In this section we determine analytically the behaviour 
of the function A{t), Eq. ([57]) . for short and long elapsed 
times t. We focus on the linear chain, when the param- 
eters are such that i/t is close to the value v^"^^ of the 
transition to the zigzag configuration. 

Short time behaviour. For short times, such that 
vtt <C 1, we expand A{t) at lowest order in the small pa- 
rameter u)y{k)t <C 1, obtaining Eq. (|40p . where F is given 
by expression (|^T|) . In the linear chain (A > 0), using 
the definition of aua in Eq. (j28[) with the corresponding 
values of the elements of the matrix TZ in Eqs. (IAip - (IA4p . 
we find 



k.(T 



4m 



N ^ 

k.cr 



UJy{k) 



(Bl) 



(B2) 



The value of F is hence proportional to the mean value 
of the transverse frequencies. 

An explicit form of the mean value can be obtained 



close to the phase transition, when A = i/t 



and for a sufficiently large number of ions. In this limit, 
we can approximate the dispersion relation Eq. (jlSp for 
small momenta q — ir/a — k as: 



A^, 



2m 



-y— 



(B6) 



The value close to the critical point and for large chains 
can be found following the same procedure applied for 
obtaining Eq. (|B4p . We get 



Ar_ 



hk\ a 
2m TT 

2'Kh 



dq 



yp+7iV 

In A -f c , 



(B7) 



where c is a constant which depends on gmax and we 
omitted terms which vanish when ^ — > 0. We now turn 
to the expression (|46|) for B{t) and with the aid of formula 
(IB3p we can approximate it as: 



B{t) = 



2h 



cos i^y{q)t 
U}y{q) 

NoiSt) , 



dq 



(B8) 
(B9) 



where we let (/max ^ oo. The function Nq{x) is the Bessel 
function of the second kind (also von Neumann function) 
[46,] . At the asymptotics, for a; ^ oo, it behaves as |46| 



Noix) 



sm X — cos X 



TTX 



(BIO) 



UJy{q) ~ \fWT}^ 



(B3) 



and hence vanishes as \l\fx 
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